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ABSTRACT 

.<^. '. 

Context. In massive-star binary systems, the interaction of the strong stellar winds results in a wind collision region between the 
stars, which is limited by two shock fronts. Besides the nonthermal emission resulting from the shock acceleration, these shocks emit 
' thermal (free-free) radiation detectable at radio frequencies that increase the expected emission from the stellar winds. Observations 

and theoretical studies of these sources show that the shocked gas is an important, but not dominant, contributor to the total emission 
in wide binary systems, while it plays a very substantial role in close binaries. 

Aims. The interaction of two isotropic stellar winds is studied in order to calculate the free-free emission from the wind collision 
' region. The effects of the binary separation and the wind momentum ratio on the emission from the wind-wind interaction region are 

a I investigated. 
Methods. We developed a semi-analytical model for calculating the thermal emission from colliding stellar winds. Assuming radiative 
(— I . shocks for the compressed layer, which are expected in close binaries, we obtained the emission measure of the thin shell. Then, we 

^ ^ computed the total optical depth along each line of sight to obtain the emission from the whole configuration. 

I \ Results. Here, we present predictions of the free-free emission at radio frequencies from analytic, radiative shock models in colliding 

O ■ wind binaries. It is shown that the emission from the wind collision region mainly arises from the optically thick region of the 

compressed layer and scales as ~ Z)"*'^, where D is the binary separation. The predicted flux density 5^ from the wind collision 
C/3 ' region becomes more important as the frequency v increases, showing higher spectral indices than the expected 0.6 value (Sy oc v", 

' where a = 0.6) from the unshocked winds. We also investigate the emission from short-period WR+0 systems calculated with our 

analytic formulation. In particular, we apply the model to the binary systems WR 98 and WR 113 and compare our results with the 
observations. Our theoretical results are in good agreement with the observed thermal spectra from these sources. 

^ . Key words. Radio continuum: ISM - Binaries: close - Stars: winds, outflows 

oo 

' 1 . Introduction In binary systems, the stellar winds of the components must 

collide, resulting in the formation of a two-shock wave struc- 

„ ■ , n , ■ .iTx.rT^ ,,T7T^s tuTe fac twcen tKe s taTS ( See , f Or Ins t ance, ElcHler (fe Uso V 1 99 3 ) . 

O Stellar wmds from hot massive stars, OB and Wolf-Rayet (WR) ^^j^^ ^-^^ ^^^^-^-^^ -^^ ^^^^^ ^^-^^ ^^^^ ^^^^^^^ 

^ ; type stars, emit free-free thermal emission detectable at ra- ^^^^.^^^ ^^^.^^.^^ ^ p.^^^^.^ 2006). The thermal compo- 

T^H . dio frequencies. Stars with spherically symmetric, isothermal, ^^^^ ^^^^^ explained as free-free emission, while the non- 

> . and stationary outflows are predicted to produce radio spectra ^j^^^^^j ^^mponent is thought to be synchrotron radiation aris- 

with a chai-actenstic frequency dependence of S . oc y" ^ (see . ^^^^ ^j^^^^^^^ accelerated at the shocks bounding the WCR 

Panagia & Felli, 1975; Wright & Barlow 1975). The specti-al wind-wind interaction zone (see, for instance, Williams et 

H index a = 0.6 results from the radial dependence of the elec- ^^^q^ ^^^^^ ^.^^ ^.^^^^^^ observations at radio frequen- 

y. tron density, n oc r -. These authors also show that variations ^.^^ ^^^^ ^ ^^^^.^j .^^.^^^ ^^^^ ^^^^ 

in this electron density behavior results in free-free spectra with ^^^.^^ nonthermal emission dominates the total spectrum. On 

a frequency dependence different from the 0.6 value. Likewise, ^^^^^ ^j^^^ ^ ^^^^^^ ^j^^ synchrotron radiation must 

Leiflierer & Robert (1991) discuss several mechanisms that may p^duced within the optically thick region of the unshocked 

cause a deviation from a 0.6 spectral mdex focusing on changes ^^^^^ ^.^^^^ ^^^^ ^^ ^^^^^ ^.^^ attenuated by 

in the ionization structure and velocity gradients in the wind. absorption (e.g.. Chapman et al. 1999; Monnier et al. 

Such effects may be present and produce small variations in the 2002). The first quantitative investigation of the effect of bina- 

obserred slope of the radio spectrum (0.6< a < 0.7). In addi- ^ity on thermalradio emission was performed by Stevens (1995), 

tion Gonzalez & Canto (20(^)8) showed that variabilities m the ^j^^ ^^^^^ ^^^^ ^.^^^^^^ ^ companion with a sti'ong wind 

wind parameters (such as velocity and mass loss rate) at injec- ^^^^^^^^^ ^^e expected thermal radio emission as compared to 

tion produce internal shocks that generate thermal continuum ra- ^ ^.^ ^^ ^^^^ ^.^^ ^^^^ ^.^^ parameters as the primary. In 

diation detectable at radio wavelengths. Their predicted spectral ^.^^ ^ ^^^^^^ ^^^^^^^ important but not dominant 

indices cleariy deviate from the expected 0.6 value of the stan- .^^tributor, while in close systems it plays a very substantial 
dard model. 
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role in the excess radio emission. In addition, Kenny & Taylor 
(2005) present colliding wind models for symbiotic star sys- 
tems by assuming good mixing of the shocked material from 
both winds. They performed three-dimensional simulations of 
bremsstrahlung radio images, and adopted fully ionized collid- 
ing winds. They found spectral differences associated with the 
viewing angle, which imply that the flux and spectral index vary 
with orbital phase. 

Pittard et al. (2006) carried out numerical models for com- 
puting the radio continuum contribution from adiabatic shocks 
in a colliding wind binary, and showed that the WCR clearly 
impacts on the thermal spectrum. They also found that the hot 
gas within the WCR remains optically thin, with an intrinsic 
emission of spectral index awcR ~ -0.1. These authors inves- 
tigated how the thermal flux from the WCR varies with binary 
separation, and found that its free-free emission scales as D^, 
where D is the binary separation. Consequently, they pointed out 
that a composite-like spectrum (which may suggest nonthermal 
emission) can result entirely from thermal processes. More re- 
cently, Pittard (2010) developed 3D hydrodynamical models for 
computing the thermal radio to submillimeter emission from ra- 
diatively colliding wind binaries (in O+O star systems; see also 
Pittard 2009). In these models, the flux density and the spec- 
trum as a function of orbital phase and orientation to the ob- 
server are investigated. They computed flux variations with or- 
bital phase, which are caused by changes in the relative position 
of the emitting components with the orbital motion. In particular, 
they investigated an eccentric system where the physical proper- 
ties (such as density and temperature) of the WCR change along 
the orbit, which is optically thick (highly radiative) at periastron 
but optically thin (adiabatic) at apastron. 

Observations of stellar winds at radio wavelengths from 
some hot stars (see, for instance, Leitherer & Robert 1991; 
Altenhoff, Thum, & Wendker 1994; Nugis, Crowther, & WiUis 
1998) show flux densities and spectral indices that differ from 
the expected value for the uniformly expanding wind model. 
On the other hand, Montes et al. (2009) presented multi- 
frequency radio observations of a sample of several WR stars, 
and discussed the possible scenarios to explain the nature 
of their emission. From these data, they found evidence for 
sources with thermal (free-free thermal emission), nonther- 
mal (dominant synchrotron emission), and composite (ther- 
mal -i-nonthermal) spectra. Nevertheless, close binaries classified 
as thermal -Hnonthermal sources cannot be completely ruled out 
as thermal sources (as pointed out by Pittard et al. 2006), and 
theoretical models are required to determine whether a compos- 
ite spectrum can be reproduced entirely by thermal emission. 

Canto, Raga, & Wilkin (1996) developed a formalism based 
on linear and angular momentum conservation for solving steady 
thin-shell problems, which is applicable to the interaction of 
non-accelerated flows. In particular, these authors found ana- 
lytic solutions to the case of two colliding isotropic stellar winds. 
Assuming that the postshock fluid is well mixed across the con- 
tact discontinuity, they obtained the shape of the thin shell, its 
mass surface density, and the velocity along the layer Canto, 
Raga, & Gonzalez (2005) also present a model for calculating 
the emission measure from thin-shell flow problems. In particu- 
lar, these authors applied the results obtained by Wilkin (1996; 
2000) and Canto, Raga, & Wilkin (1996) for a stellar wind in- 
teracting with an impinging ambient flow with a density strat- 
ification. From their model, simple predictions of the free-free 
emission can be made. 

In this work, we study the case of the interaction of two 
spherically symmetric stellar winds. We present a formulation 



for calculating the free-free emission from the wind- wind inter- 
action region, with which the total thermal spectrum of colliding 
wind binaries can be obtained. This paper is organized as fol- 
lows. In § 2, we describe the model. In § 3, we show the results of 
our model for the predicted thermal fluxes and the corresponding 
spectral indices for massive binary systems. A comparison with 
observations of this kind of sources are presented in § 4. Finally, 
we give our conclusions in § 5. 

2. The analytical model 

We consider a binary system that contains colliding isotropic 
stellar outflows. This situation is shown in a schematic way in 
Figure 1 . We assume two ionized winds that move radially away 
from the stars separated by a distance D. The stellar winds are 
generally accelerated to hypersonic speeds and reach a good 
fraction of their terminal velocities in a few stellar radii. In 
our model, we assume that the winds are ejected with termi- 
nal speeds. The interaction of these two outflows gives rise to 
the wind collision region (WCR) that is, in principle, composed 
of two shocks separated by a contact discontinuity (Luo et al. 
1990). If the shocks are radiative (post-shock cooling becomes 
very efficient leading to a large compression), the shock fronts 
collapse onto a thin shell and the width of the WCR can be ne- 
glected. 

Let nil and vi be the mass loss rate and the velocity of the 
wind source at the origin of the spherical coordinate system 
(R, 0, (p), and m2 and V2 the corresponding quantities for the wind 
from the source located at a distance D. Using the geometric re- 
lation 

R{6) ^ DsinOi CSC (6 + 61), (1) 

where R{6) is the radius of the layer. Canto, Raga, & Wilkin 
(1996) obtained the exact (implicit) analytic solution 

6*1 cot 6»i = 1 -1-^8(6' ctg 61- 1), (2) 

where the dimensionless parameter j3 = (mi vi)/(m2 V2) is the 
wind momentum ratio. For a given 6, equation (|2| can be solved 
numerically for 0i , and the radius is then obtained from equation 
([T)- For the case of an isotropic stellar wind, the flow is axisym- 
metric, and therefore the spherical radius R does not depend on 
the azimuthal angle 0. In addition, the stagnation point radius, 
Rq = j6'^^D/(l H-jS'^^) is obtained from the ram-pressure balance 
condition (Canto et al. 1996). 

The pressure within the shell just behind each shock can be 
written in terms of the preshock density, and the velocity com- 
ponent normal to the shell. The normal vector to the surface 
of the shell can be obtained from the gradient of the function 
F - r' - R{0), where r' is the radial coordinate, and the bow 
shock is described hy F -Q. In this way, h - VF/ | VF |, where 
V denotes the vector differential operator gradient in spherical 
coordinates (see Wilkin 1997). For isotropic outflows dR/d4> - 
(see eq.([T]i), and then 

^ R-l/R(dR/d6)6 
V(l + l/R^ {dRIdOf) 

where R and 9 are unit vectors measured along the directions R 
and 9, respectively. 

On the other hand, the velocity vectors of the stellar winds 
are given by 

vi = vi^, (4) 
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Fig. 1. Schematic diagram showing the interaction of two spher- 
ical winds that move radially away from the stars. Source 1 of 
the wind with velocity vi is located at the origin of the coordinate 
system, and source 2 of the wind with velocity vo at a distance 
D along the symmetry axis. The shape of the layer where the 
winds collide (at Ro along the symmetry axis) is given by the 
curve R{d). The angles 9 and 9\ are measured from the positions 
of the stars to the intersection point between a line of sight (with 
impact parameter r) and the layer. We assume azimuthal sym- 
metry (no dependence on angle (p). The observer is located in 
the orbital plane along the z-axis of symmetry. 



and 

V2 = -V2 COS {0 + 9\)R + V2 sm(9 + 6i)0. 



(5) 



From equations ([3]l, © and (|5]l it follows that the normal 
component at every point (R, 6, (p) of the pre-shock velocity of 
the stellar winds are obtained by 



R 



V\,n = Vl 



V2,n = -V2 ■ 



yjR^ + idR/d9)^ 

R cos i9 + 9i) + {OR/ 89) sin (9 + 9i) 



^/WT{^RJ^^ 



(6) 



(7) 



As a result, the pressure within the thin shell just behind each 
shock front can be written as 



^1 =Pi,ovto/.v,i(0,0i), 
and 



(8) 



(9) 



where we have used the wind parameters pi o and Vj q at the stag- 
nation point Ro to nondimensionalize the equations, and 

R'- 



and 



Pi 



(10) 



/.,2(0,0l) = (-^^|X 
\P1,0 V\ Qj 

{R cos {9+9i) + idR/d9) sin {9 + 9i)f 
R^ + {dRld9f 

where p\ and p2 are the preshock densities of the stellar winds 



(11) 



2. 1. The emission measure of fhe fh/n she// 

Let / be a longitude-coordinate measured inwards from the wind 
shock of source 2 and perpendicular to the thin shell (as shown 
in Fig. 2). Then, the shock fronts are located at Z = and I = h, 
where h is the position-dependent thickness of the shell. We as- 
sume that the post-shock gas is well mixed, so that it flows along 
the shell at a velocity v, which is independent of I. If the whole 
flow is photoionized and approximated as an isothermal flow, the 
pressure as a function of the parameter I can be calculated from 
the hydrostatic equation 



dP 
ll 



-P8^ 



(12) 



where P = pc^ (with c, being the isothermal sound speed), and 
the centrifugal acceleration g = v^/Rc, with Rc the radius of 
curvature of the thin shell. Integrating equation (IT2t with the 
boundary condition f 2 = P{1 - 0), we obtain 



P(0=p(0< = P2e 



-l/H 



(13) 



where H - c^Jg is the scale height parameter of pressure (or 
density). It follows from equation (fT3l l that the surface density is 
given by 



cr 



-f 

Jo 



p(l)dl^^(P2-Pl). 



(14) 



Using equations ( fT3b and (fT4l) . we can now calculate the emis- 
sion measure of the thin shell. 



EM 



-f 

Jo 



(p/m) dl = 



cr 



■(P2 + P1), 



(15) 



where m - finip is the average mass per particle of the photoion- 
ized gas, /i the mean atomic weight per electron, and irip the mass 
of the proton. 

Equation ( fTSl l gives the emission measure of the thin shell 
in terms of its surface density. Next, we must obtain the sur- 
face density as a function of position on the shell. We first note 
that the shocked gas can be described as a flow of streamlines 
with constant azimuthal angle <p. The mixing of the two shocked 
winds is assumed to be instantaneous, so that the flow will have 
a unique velocity v at any location with in the shell. The mass 
(per unit time) AM flowing along a ring of radius r of the layer 
is then given by 



AM = 2KRsm9o-v, 



(16) 



which must be equal to the mass injection rate (into the same 
solid angle) of the stellar winds, that is. 



AM= ^/„,(0,0i), 
with 



f,ni9, 0i) = 1 - COS + -^(1 - COS 0i) . 

j6v2 



(17) 



(18) 



Based on considerations of linear momentum conservation, 
it can be shown that the velocity along the shell (see also Canto, 
Raga, & Wilkin 1996; Canto, Raga, & Gonzalez 2005) can be 
written as 



V = Vl - 



^fi(9,9i)+fi(9,9i) 



(19) 
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Fig. 2. Schematic diagram showing a thin shell (limited by two 
shock fronts) resulting from the interaction of two stellar winds. 
The shocked gas flows along the shell at a velocity v and has a 
surface density cr. The coordinate I is measured inwards from 
the wind shock of source 2, normal to the locus of the thin shell. 



where the functions /,. {9, 9\) and /, {6, 6i) are deduced from r- 
momentum and z-momentum (being r and z the directions of the 
cylindrical radius and the symmetry axis, respectively), and are 
given by 



1 



me,) 



and 



9 - sin ffcos 9 + —{9i - sin 9\ cos 9i) 



sm 



■ - sin^ 6»i 



(20) 



(21) 



Substitution of equations (fT7li- (l2TI) into equation ( fTSI l gives 
the surface density, 

CT^(rQU9,9i), (22) 
where ctq = mil{2nfiDv\) and 

f^(9, 6»i) = ^ sin (6 + 9i) esc 6i esc 9 ' ^' .(23) 

^f^(e,90+f?{9,9:) 

Finally, it follows from equation (flST l that the emission mea- 
sure of the thin shell (as a function of position) is given by 



EM(9,9i) = EMoU(9,ei) (9,61) + A,2 (9,9i)] , 
with EMo = (Topifl v\j{2np-cl). 



(24) 



2.2. Predicted thermal radio emission from a binary system 

To calculate the radio-continuum emission from the whole sys- 
tem, it is necessary to compute the total optical depth along each 
line of sight, which will have the contribution from the stellar 
winds and also from the thin shell. Then we estimate the in- 
tensity emerging from each direction and calculate the flux by 
integrating the intensity over the solid angle. 

Given the emission measure EM(9, Ox ) of the shocked layer 
(eq.(l24l)). the optical depth perpendicular to the WCR is calcu- 
lated by TwcR.±ifi,9i) = EM{9,9i)x{v), where xiv) =8.436 x 
10"^ v"^ ' with the frequency v in Hz (for an electron tempera- 
ture Te - lO'^K). By defining a critical frequency, such that 



-2.1 



i,47rmv-^y 



it can be shown from equation (l24l l that 



j/r(0, Oy) [/„..i (0, Oi) + U,2 (9, 9i)] ., 



(26) 



where D = D/Rq. For a line of sight intersecting the thin shell at 
an angle cos"^ (z-h) from the normal (being z = cos 9R-sm99), 
the optical depth is then given by 



Rcos9 + sm9(dR/d9) 



(27) 



where R is in units of Rq. 

Let us now calculate the contribution of the unshocked stel- 
lar winds to the optical depth. Consider a line of sight that in- 
tersects the thin shell at a point (rctgfl, r), as shown in Figure 
1. According to Panagia & Felli (1975) and Wright & Barlow 
(1975), the optical depth along the line of sight of the wind 
source located at z = is obtained by 



T.,l(9) 

where 



r-rCtgi 

%J —00 



n„A{z)xiv)dz, 



(28) 



(29) 



where nio - pifl/fh (- mi/4nmviflRQ ) is the number density of 
the flow at the stagnation point. From equations (l25T l. (l28T l, and 
(|29] l, it follows that 



with 



-2.1 



(30) 



ag9 



n 

2 ctg^ 9 + 1 



Analogously, the optical depth of the wind source located at z = 
D is calculated by 



T„,2(0)^ nl^{z)xiy)dz, 

JrCtge 

where 



«H.,2(2) = n2,{) — 



1 



[z - 6]2 + r2 ' 



(31) 



(32) 



where n2,o = Pifi/in (= nii /4-mnv2fiRl) is the number density of 
the wind at the stagnation point. Substitution of ( l32b into equa- 
tion OU gives 



-2.1 



/v,a(0) , 



(33) 



with 



(25) /„,2(0) 



OT2 Vl 



m\V2j 2r^ 



1 



{fcig9-D)f 
+ (rctg 61-5)2 
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10^ km s-\ mi = 1.25 x 10"^ Mq yr 



0.010 



0.100 



Fig. 3. Predicted free-free emission from a colliding wind binary. 
We have adopted the parameters vi = 10^ km s"', nii = 1.25 x 
10-^ Mo yr-i, and = 10^ km s"', m = 5 x lO^^ Mg yr"' 
for the wind sources. The binary separation is set to D = 4 AU. 
The dotted and dashed lines represent the fluxes (oc v*'-^) from the 
wind sources 1 and 2 (see also Fig. 1.), respectively. The intrinsic 
thermal emission from the WCR (dot-dashed line) and the total 
emission (solid line) from the binary system are also shown. The 
behavior of the curves is described in the main text. 



- arctg 



'fctgd- D 



We assumed that the system is far enough (D <K L, being L 
the distance to the observer) so that the lines of sight intersecting 
the central stars can be ignored. Then, the radio-continuum flux 
from the colliding wind binary can be calculated by 



Sy = 2nBy 



[l-e 



-T(e,ei) 



] fdr. 



(34) 



where t(0,0i) - TwcRiO,0\) + ^K^xiff) + T„;i{ff) is the total op- 
tical depth along the line of sight, By{- IkTv^jc^ with k the 
Boltzmann's constant, and c the light speed) is the Planck func- 
tion in the Rayleigh-Jeans approximation, and f{9oc) is the im- 
pact parameter at the asymptotic angle 0oo (corresponding to 
7? — > oo) of the thin shefl, which can be found from 

TT 

6oo ~ tan 



1-yS 



(see Canto, Raga & Wilkin 1996). Defining the parameter - 
4nkT(R()/ cL)^v^., we finally obtain 



4- 



Jo 



[l-e 



-T(m h 52 



R'-(e)smecos6de- 



f 

Jo 



[1 - e"^(*>] ^(6*) sm^6(dR/d6) dO 



(35) 



Equation (l35T l represents an self-similar solution for the free-free 
radiation from colliding stellar winds from massive binary sys- 
tems. 

We apply the model for computing the emission from a col- 
liding wind binary using stellar wind parameters similar to those 
from a typical WR-hO binary system. That is, we have adopted 



for the O type star, 
1 - 5 x 10"^ Mp 



and for the WR component, V2 = 10 km s ', m2 - j lu iviq 
yr"'. In addition, we have assumed a separation D =4 AU be- 
tween the stars (so that the WCR is radiative; see § 3). The result 
is shown in Figure [3] It shows the contribution to the total radio 
emission of the unshocked stellar winds and the WCR. As ex- 
pected for ionized stellar envelopes (see § 1), the flux densities 
from the stellar winds increases as ~ v'"'. However, deviation 
from the expected 0.6 value of the spectral index at high fre- 
quencies is observed from the emission of the wind source 1. 
This is probably due to the presence of the WCR inside the op- 
tically thick region of the wind. We note that as v increases the 
emission from the WCR becomes more important. The intrinsic 
flux density from the WCR shows a spectral index of ~ 1.1, con- 
sistent with the numerical models (in massive O+O type binary 
stars) developed by Pittard (2010). At low frequencies, the total 
flux density grows as 5 y ~ v"'^ approaching the emission from 
the stronger wind. On the other hand, at higher frequencies, the 
radio spectrum from the system approaches the flux density from 



the WCR (5v 



,1.1 



3. Thermal radio-continuum emission from 
colliding wind binaries 

As mentioned in § 1, radio observations (e.g. Moran et al. 1989; 
Dougherty, Williams & Pollaco 2000; Montes et al. 2009) and 
theoretical models (see, for instance, Dougherty et al. 2003; 
Pittard et al. 2006; Pittard 2000) of massive binary stars have re- 
vealed strong shocks formed in the wind-wind interaction zone. 
These shocks emit radio-continuum radiation consistent with 
thermal emission that may produce variations in the expected 
flux densities and spectral indices from ionized stellar winds. In 
this section, we apply the model of colliding-wind binary sys- 
tems developed in § 2 by adopting different wind parameters 
of the components. Our model assumes strong radiative shocks 
which collapse onto a thin shell, so that the width of the WCR 
can be neglected. We investigate the contribution of the WCR to 
the thermal radiation from colliding-wind binaries, and the effect 
of the binary separation on the radio spectrum. 

Stevens, Blondin & Pollock (1992) investigated the collision 
of stellar winds in early-type binary systems. They studied the 
role of radiative cooling in the structure and dynamics of col- 
liding wind binaries. In this work, the importance of cooling 
in a particular system is quantified using the cooling parame- 
ter (» 0.15 VjC/au/^s, where V3 the wind velocity in units of 
10^ km s"' , dp,\j the distance to the contact discontinuity in units 
of AU, and M5 is the mass loss rate in units of 10"^ Moyr"'), 
defined as the ratio of the cooling time of the shocked gas to the 
escape time from the intershock region. For models with x ^ 
I, the postshock flow can be assumed to be adiabatic, while it 
is roughly isothermal for models with x ^ ^ - Numerical simu- 
lations by these authors show that the cooling result in the for- 
mation of a thin dense shell, confined by isothermal shocks, as 
X drops below unity. Antokhin, Owocki, & Brown (2004) have 
also investigated the narrowness of the cooling layer using the 
ratio IqIR (where /q is the cooling length and R the radius from 
the star source), which is closely related to the above parameter 
X- These authors show that the ratio IqIR can serve as well to dis- 
tinguish between adiabatic and radiative shocks, with /q//? > 1 
implying an adiabatic shock and l^jR < 1 a radiative one. 

Nonetheless, radiative pressure that moderates the wind- 
wind collision may be important in close hot-star binaries. An 
initial analysis was developed by Stevens & Pollock (1994), who 
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Fig. 4. Optical depth TwcRi^^ GHz) of the WCR as a function of impact parameter r (left panels), and the predicted flux density 
Sy at radio frequencies (right panels) for the models of Table 1 . In top panels, the dot-dashed, solid and dotted lines represent the 
optical depth and the radio spectrum from the models Bl, B2, and B3, respectively. In bottom panels, the corresponding values of 
models B2, B4, and B5 are shown by the solid, dot-dashed, and dotted lines, respectively. The physical description of the plots are 
given in the text. 



investigated the dynamics of colliding winds in massive binary 
systems. They show that, in close binaries, the radiation of a lu- 
minous star inhibits the initial acceleration of the companion's 
wind towards the stagnation point. These radiative forces re- 
sult in lower velocities than those expected in single-star models 
moderating the wind collision. In other work, Gayley, Owocki 
& Cranmer (1997) investigated the potential role of the radia- 
tive braking effect, whereby the primary wind is decelerated by 
radiation pressure as it approaches the surface of the compan- 
ion star. These authors conclud that radiative braking must have 
a significant effect for wind-wind collision in WR+O binaries 
with medium separations {D < 0.5 AU). Furthermore, Parkin 
and Pittard (2008) carried out 3D hydrodynamical simulations 
of colliding winds in binary systems. They shows that the shape 
of the WCR is deformed by Coriolis forces into spiral structures 
by the motion of the stars. The shape of the shock layer is more 
deformed in systems with eccentric orbits. In addition, Pittard 
(2010) developed 3D hydrodynamical models of massive O-i-O 
star binaries for computing the thermal radio to submillimeter 
emission. In this work, flux and spectral index variations with 
orbital phase and orientation of the observer (from radiative and 
adiabatic systems) are investigated. These authors found strong 



variations in eccentric systems caused by dramatic changes to 
the physical properties (density and temperature) of the WCR, 
which is radiative (optically thick) at periastron, but adiabatic 
(optically thin) at apastron. 

3.1. The predicted radio emission from radiative sliocks in 
colliding wind binaries 

Here, we present analytic predictions of the thermal spectra at 
radio frequencies from different radiative models, which satisfy 
the condition of radiative shocks (x < I) for the WCR. We have 
estimated the value of ^ for each shock of the WCR. In Table 1, 
we list the different scenarios of colliding wind binaries that we 
have studied in this paper In models B1-B3, we have assumed 
identical wind sources (J3- 1) in order to investigate the effect of 
the binary separation. On the other hand, in models B2, B4, and 
B5, we have assumed different wind momentum ratios for the 
binary systems, while the distance between the components is 
unchanged. In models B4 and B5, we indicate the highest value 
of^ for each system. These models for the radio emission from 
colliding wind binaries are presented in Figure We show the 
optical depth of the WCR as functions of the impact parameter 
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TABLE 1 

Parameters of the colliding winds 



Model 


P 


D/AU 


X 


RalD 


Bl 


1 


8 


0.12 


0.5 


B2 


1 


4 


0.06 


0.5 


B3 


1 


1 


0.02 


0.5 


Model <** 


P 


O/AU 


X 


RolD 


B4 


0.25 


4 


0.15 


0.33 


B5 


0.025 


4 


0.65 


0.14 



(a) vi,2 = 10' km s"'; m,,2 = 5 x 10"^ Mq yr"' 

(b) vi'z = 10' km s-'; mz' = 5 x 10^^ Mq yr"'; m, = 1.25 x 10^^ M, 
yr"' for model B4 and mi = 1.25 x 10"'' Mq yr"' for model B5. 



r (left panels) and the predicted flux density at radio frequencies 
(right panels). Top panels show the results of our Models BI- 
BS, while the models B2, B4, and B5 are shown in the bottom 
panels. All models are calculated for an observer located in the 
orbital plane along the symmetry axis (which corresponds to a 
system with an inclination angle / = 90° with the most powerful 
stellar wind in front, as shown in Figure 1). 

First, we observe that, for low-impact parameters (r//?o ^ 
1), the optical depth tvtcs (43 GHz) of the WCR (Figure |4} 
a) does not depend on the value of r and scales as D^^ . On 
the other hand, for higher impact parameters (r/7?o s> 1), 
^1^^(43 GHz) oc r^^ in all models. At a given line of sight, we 
also note that, in this limit, the optical depth increases as D^. The 
transition impact parameter is r oc Rq (see, also. Appendix A). 
This behavior of the optical depth (and therefore of the emission 
measure of the thin shell) can be explained as follows (see Figs. 
IA.ll and [A.2l l. We consider a given impact parameter r and vary 
the binary separation. For systems with r/7?o » 1 (t, oc r^^), as 
the binary separation (and therefore Rq) increases, the ratio rjRo 
decreases. The normal components of the shock velocities (for 
the two shock fronts that bound the WCR) increase for higher 
values of Rq, which result in higher pressures within the shocked 
layer (eqs. ||8] and ||9l). Thus, both the emission measure (eq. 
ifTSl ) and the optical depth (eq. Il27l ) of the thin shell increase. 
On the other hand, for systems with r//?o 1 {ts + '!".s['"/^o]), 
as R<^ increases the preshock densities (for the two shock fronts 
of the WCR) decreases, which result in lower pressures within 
the shocked layer Consequently, the emission measure and the 
optical depth of the WCR diminish. 

Second, we show the predicted flux density at radio frequen- 
cies for models B1-B3 (Figure|4]-b). We note that, at frequencies 
V > 1 GHz, as the binary separation increases, the emission from 
the binary becomes more important. In order to explain this be- 
havior of the spectrum, we investigated the emission from the 
shocked layer as a function of the impact parameter r (see §A.3 
of the Appendix A). Our results are presented in Fig. IA.3I We 
show that the main contribution to the flux comes from the opti- 
cally thick region of the shell, so the emission from the optically 
thin region can be neglected. In addition, our model predicts that 
the emission from an optically thick WCR (which is bounded by 
radiative shocks) scales as S'y oc D'*'^ (see Fig. lA.4l ). This result 
contrasts with the inverse dependence with the binary separation 
described by Pittard et al. (2006) for an adiabatic and optically 
thin WCR, where the emission is expected to scale as D ' . 

Finally, as mentioned above, the effect of varying the wind 
momentum ratio yS on the value of TiycR(43 GHz), and the total 
spectrum, is investigated from models B2, B4, and B5 (bottom 
panels of Fig. |4]i. As we expected from equations (|26]|-(|27T|. it 
can be observed that the optical depth of the WCR increases 



TABLE 2 

Wind parameters of colliding wind models 



Model 


/3 


D/A(/ 


X 






WR 98 


0.20 


0.5 


0.6 


0.31 


0.6 


WR 113 


0.21 


0.5 


0.6 


0.31 


0.6 



with the mass loss rate m\ of the secondary component (Fig. 
|4]-c). In addition. Figure |4]-d shows the predicted thermal radio 
spectra from these models. The plots show how sensitive our 
analytic model is to the wind momentum radio, fi. We note that, 
at high frequencies, the total flux density S y increases with the 
parameter p. 

4. Comparison with observations 

As we mention in previous sections, the model presented in this 
paper assumes a thin shell approximation for the colliding wind 
region. In the presence of efficient cooling (^ < 1), the WCR is 
confined by radiative shocks, which are expected in close bina- 
ries. Consequently, we investigate the emission from short pe- 
riod (< 1 yr) WR+O systems calculated with the analytic model 
developed in § 2. Besides, radiative breaking can be important 
for binaries with a medium separation D < 1007?© (orbital peri- 
ods ~15 days), which prevents the applicability of our model for 
systems with shorter periods. 

Recently, Montes et al. (2009) have reported observations of 
the flux density and spectral indices for the close systems WR 98 
113, 138, and 141. In this section, we apply the analytical model 
to the binaries WR 98 and WR 113 and compare them with ob- 
servations. Table 2 lists the adopted stellar parameters of these 
sources in our colliding wind model, which suggest the presence 
of radiative shocks in the wind-wind interaction zone. We indi- 
cate the highest value of x for each binary system. We indicate 
the highest value of x for each binary system. We show that the 
contribution of the WCR to the total emission gives a possible 
explanation to the observations. 

WR 98. This close binary has been identified as a double line 
binary (WN7o/WC+08-9) with an orbital period, P = 47.8 days 
(Gamen & Niemela 2002). Abbott et al. (1986) classified WR 98 
as a nonthermal radio source. Recently, multi-wavelength obser- 
vations (from 5 to 23 GHz) of this source have revealed that the 
spectral index changed from ~ 0.26 to ~ 0.64 in a period of time 
of ~ 15 days (Montes et al. 2009). This behavior was interpreted 
as indicating a binary influence over the radio spectrum, pos- 
sibly from a variable nonthermal contribution that is absorbed 
at certain orbital phases, turning the spectrum into a "thermal 
state". The stellar parameters of this binary indicate that a radia- 
tive WCR (x ~ 0.6) is formed within this system. 

We applied the model to the binary system WR 98 to inves- 
tigate if a thermal component of emission from the WCR is able 
to contribute significantly to the total spectrum. Our results from 
the model are compared with the observations at the thermal 
state (a ~ 0.64) reported by Montes et al. (2009). We assumed 
for the WR star a mass loss rate M-wr = 3.0 x 10"^ Mq yr"' (up- 
per limit derived from radio observations at the thermal phase 
at 8.4 GHz when a ~ 0.64; Montes et al. 2009), and an ejec- 
tion velocity vwr = 1200km s"' (Eenens & Williams, 1994). 
The best fit to the observations was found from the O star pa- 
rameters Mo - 4.0 x lO^^Moyr"' and vq - 1800 km s"'. We 
have assumed a mean atomic weight per electron fi = 4.2 and 
an average ionic charge Z = 1 . 1 for the WR wind, and fi - 1.5 
and Z - I for the O star (which are typical values for O-type 
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1 00.00 r 




0.10 r 



0.01 V'-'". . . 

1 10 100 

V (GHz) 

Fig. 5. Comparison between our model and the flux densities at 
5, 8.4, 23 and 250 GHz ('+' sign; uncertainty in the flux mea- 
surements are given by the length of the bars) of the binary 
system WR 98 obtained by Montes et al. (2009) and Altenhoff" 
et al. (1991), respectively. We have assumed the stellar wind 
parameters Mvi/j? = 3.0 x lO^^Moyr"', v^^r - 1200 kms"'. 
Mo = 4.0 X 10~''Moyr-', and vq = 1800 km s '. The stai's ai-e 
separated by a distance D ~ 0.5 AU. The dotted and dashed 
lines represent the fluxes (oc v°^) from the winds of the WR and 
O type stars, respectively. The radiation from the WCR (dot- 
dashed line) and the total flux density (solid line) from the bi- 
nary source are also shown. The physical description of the plot 
is given in the text. 



100.00 r 




0.01 L . ^^^^^ . ^ . J 

1 10 100 

V (GHz) 

Fig. 6. Comparison between our model and the flux densities at 
5, 8.4, 23 GHz ('-n' sign; uncertainty in the flux measurements 
are given by the length of the bars) of the binary system WR 113 
obtained by Montes et al. (2009). We have assumed the stellar 
wind parameters MvvR = 2.0 xlO"^ Mq yr"', wt's - 1700km s"'. 
Mo = 0.4 X lO^^Moyr^', and vo = 1800 kms"'. The stars are 
separated by a distance D ~ 0.6 AU. The dotted and dashed 
lines represent the fluxes (oc v*'^) from the winds of the WR and 
O type stars, respectively. The radiation from the WCR (dot- 
dashed line) and the total flux density (solid line) from the binary 
source are also shown. The physical description of the plot is 
given in the text. 



stars; see Bieging et al. 1989). Using equation (19) in Canto et al. 
(1996) and considering that most of the emission from the WCR 
arise from impact parameters, r, such that < Oca (see Table 2), 
we can estimate the ratio of mass entering into the WCR from the 
WR stellar wind, to that from the O star, AMwr/AMo > 0.75. 
Thus, we assume that the material within the shock is mainly 
composed of what comes from the strongest wind, assuming the 
same values of /i and Z for the shocked material as those used for 
the WR wind. From the value of a sin ; ~ 100 Rq, (where a is the 
semi-major axis and / is the inclination angle of the orbit) deter- 
mined by Gamen & Niemela (2002), we assume D ~ 0.5 AU for 
the binary separation, which is the maximum value for a (when 
i ~ 90°). Since this value represents an upper limit for D and 
because the geometric configuration represents the WR stellar 
wind in front of the O star, we are determining the minimum 
expected for the WCR thermal contribution. 

In Figure |5] we show the predicted spectrum from WR 98 
system. Observational data (^sghz = 0.58 +0.06mJy, 5'8.4ghz - 
1.18 + 0.05 mJy, and SiiCHz = 1-94 + 0.15 mJy) by Montes et 
al. (2009) are also plotted. It can be seen that our model pre- 
dicts an increase in the flux density at high frequencies, which 
results in steeper spectral indices than the expected 0.6 value of 
the stellar winds. Also the WCR emission becomes comparable 
to that from the WR wind at ~ 70 GHz, and at a frequency of 
250 GHz, the shock layer produces an excess of emission of a 
factor > 2 over the value expected for the WR wind (8 mJy). 
To test this prediction of the model at higher frequencies, we 
included the flux density 5250GHz = 19 + 5mJy measured by 
Altenhoff et al. (1991), which seems to agree with the 22 mJy 
predicted for the total flux density at 250 GHz. Furthermore, the 
spectral index of the WCR emission , q-wcr ~ 1 -2, is similar to 



what is derived from the observations at this frequency range, 

ff22.5-250GHz ~ 0.95. 

WR 113. This source was detected at radio frequencies 
(8.4 GHz) for first time by Cappa et al. (2004). The period of this 
spectroscopic binary (WC8-I-08-9) is P = 29.7 days (Niemela 
et al. 1996). Recently, simultaneous multifrequency observa- 
tions (from 5 to 23 GHz) of WR 113 by Montes et al. (2009) 
showed an spectral index a ~ 1.06. We have modeled this bi- 
nary adopting the wind parameters Mwr - 2.0 x lO^^Moyr ' 
(Lamontagne et al. 1996) and vwr = 1700 km s"' (vanderHucht 
et al. 2001) for the WR star, and Mo = 0.4 x lO^^Moyr"' 
and v'o = 1800 km s"' for the O type component, separated by 
a distance D - 0.6 AU (from the value a sin/ - 129/?©, and 
/ ~ 70 determined by Lamontagne et al. 1996). We assume a 
mean atomic weight per electron fi - 4.7 and an average ionic 
charge Z - l.l for the stellar wind of the WR wind (from Cappa 
et al. 2004), and ju = 1.2 and Z = 1 for the O star As in 
the case of WR 98, we assumed that the shock is mainly com- 
posed of material from the WR wind (AMwr/AMo > 0.65) , 
and the same fi and Z values were used. Figure |6] shows the 
predicted spectrum by the model, as well as observational data 
(55GHz = 0.22 + 0.03 mJy, 58.4GHz = 0.47 + 0.04 mJy, and 
523 GHz = 1.27 + 0.07 mJy) from the WR 113 system by Montes 
et al. (2009). We note from the figure that the contribution of the 
WCR starts to dominate the emission from the system at a lower 
frequency (v ~ 9 GHz) than in the case of WR 98, mainly be- 
cause the flux density from the WR star is ~ 3 times lower than 
for WR 98. At a frequency of 250 GHz, the flux density from 
the thin shell becomes a factor of ~ 6 greater that the expected 
value from the unshocked winds. Observations at such high fre- 
quencies are required in order to verify the high spectral indices 
predicted by our model. 
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From our model, we have shown that a thermal contribution 
from the WCR in WR 98 and WR 113 is Ukely to be detected 
at high frequencies. However, as we pointed out, we can only 
present a lower limit for this contribution, owing to the uncer- 
tainty of parameters such as the binary separation, D. On the 
other hand, Pittard (2010) have investigated the variability in the 
flux density (from radiative systems) due to the orbital motion 
of the stars and orientation of the observer. The configuration 
adopted in our model (with an inclination angle of 90° and or- 
bital phase ~ 0) corresponds to the minimum contribution to the 
total thermal emission expected from radiative shocks in binary 
systems with an inclination angle ~ 90°. In this way, the excess 
of emission predicted here is expected to be variable and mod- 
ulated by the orbit motion. WR 98 was reported as a variable 
source by Montes et al. 2009, changing its spectral index from 
~0.26 to ~0.64 in a period of ~15 days. Although the flat spectral 
index (~0.26) resembles the one predicted by Pittard et al. 2006 
for an adiabatic shock, this close system is likely to be radiative, 
and such behavior was explained as result of a nonthermal con- 
tribution escaping the absorption. On the other hand, for WR 113 
there are no observations at high frequencies or clear evidence of 
variability to support the contribution predicted here; therefore, 
high-frequency observations are required to confirm it. 

5. Summary and conclusions 

In this paper, we have presented an analytical model for calc ulat- 
ing the thermal contribution at radio frequencies from the WCR 
in massive binary systems. The strong stellar winds of both 
components collide and form the WCR (bounded by two shock 
fronts) between the stars, which emits both thermal (free-free) 
and nonthermal (synchrotron) radiation. In particular, the non- 
thermal component is expected to be highly absorbed in close 
binaries, so the thermal flux must be the dominant contributor. 
In addition, particle acceleration should not occur if the shocks 
are collisional such as in short period systems where the gas den- 
sity is high. 

Based on considerations of linear and angular momentum 
conservation, we developed a formulation for calculating the 
emission measure of the WCR, from which simple predictions of 
the thermal emission can be made. We assumed that the WCR 
is bounded by radiative shocks, so that the width of the WCR 
can be neglected and the wind-wind interaction zone can be ap- 
proximated as a thin shell. Since highly radiative shocks are ex- 
pected to occur in short period systems, we assumed difl'erent 
scenarios of colliding wind close binaries. Using typical parame- 
ters of massive WR-i-O binaries, we showed that the compressed 
layer emits free-free radiation that plays a very substantial role in 
the emission from the whole system. The radio-continuum spec- 
trum obtained by our model clearly deviates from the behavior 
>Sv oc v''-^ predicted by the standard stellar wind model. 

In previous models (Pittard et al. 2006), the expected ther- 
mal emission from binary systems where the WCR remains op- 
tically thin (confined by adiabatic shocks) was investigated. In 
this case, the flux density scales as D'^; consequently, the emis- 
sion from the WCR increases as the binary separation decreases. 
In this work, we study the thermal radiation from systems where 
the WCR is highly radiative. The effects of the binary separation 
and the wind momentum ratio on the total spectrum are inves- 
tigated. Our model's results show that the relative contribution 
of the WCR to the total emission mainly arises from the opti- 
cally thick region of the layer, which flux density scales as 
(where D is the stellar separation). Thus, the impact of the WCR 
to the total spectrum (in contrast to the optically thin case) be- 



comes more important as the binary separation increases (until 
the WCR becomes optically thin). 

Finally, we applied the analytical model to massive binaries, 
and calculated the thermal radio-continuum emission from the 
short-period WR+O systems, WR 98 and WR 113. Assuming 
a set of wind parameters consistent with observations of these 
sources, we find that the WCR must be confined by radiative 
shocks, and so the thin shell approximation can be applied. In 
both sources, the compressed layer generates thermal radiation 
that produces (at high frequencies) an excess of emission over 
the expected values from the stellar winds. Comparison with re- 
cent observations (by Montes et al. 2009) from these sources 
showed that our model can satisfactorily reproduce the flux den- 
sities and spectral indices. 
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Fig. A.l. Nondimension emission measure as function of rjRd. 
The dotted and dashed lines represent the functions fs of the 
shocked layer, and /„, of the stellar winds, respectively. The total 
emission measure /, + is also shown (solid line). The physical 
description of the plot is given in the text. 



Fig. A.2. Nondimension emission measure as function of r//?o- 
The function {r/Roy^ fs (dotted line) of the shocked layer, and 
(r/Roy^ fw (dashed line) of the stellar winds are presented. The 
solid line represents the total emission measure (/IRqY^ (fs+fw)- 
The physical description of the plot is given in the text. 



Wilkin, F. P. 1997, Ph.D. Thesis, 

Wright, A. E., & Barlow, M. J. 1975, MNRAS, 170, 41 



Next, we combine equations (fT6b and (IA.5b in order to find 
the surface density, 



Appendix A: Emission from a binary system with 
identical stellar winds 

We consider two stars with identical winds. The stellar winds 
are ionized and isotropically ejected with terminal velocity v (- 
vi = V2) and mass loss rate m {ihi - 1112). In this particular case, 
the wind momentum ratio of the interacting winds fi - \, the 
stagnation point radius Rq - Djl, and the angle 6-61 (see Fig. 
1). 

First, we estimate the emission measure of the shocked layer 
as follows. From equations ( fTSl l, ( l20l l, and ( 1211 1, one obtains 



/;„(0) = 2(l-cos0), 



(A.l) 



m\{l - cos fff- 



nv\Rsm9(0 - smO cos 6) ' 



(A.6) 



where we have used the velocity v given by equation (IA.4I I. 

In this simple case R - Rq/cos 9, and the normal component 
of the preshock velocity of the stellar wind vi „ - v\ cos 9 (see 
eq. ||6l). Therefore, the pressure Pi 
just behind the shock front can be written as 



Pi v? within the thin shell 



4jrRl 



(A.7) 



Note from equations O and (|9]l that Pi - P2, as expected. From 
equation (fTSl l, it follows that the emission measure of the shell is 
EMs - crPi/rfP' c^. Using equations (IA.6I) and (IA.7b , it can be 
shown that. 



f,.(9) = 0- sine cos 0, 
and, 

/.(0) = o. 



(A.2) 



(A.3) 



Substitution of eqs. ( IA.ll )-( lA3l l into eq. ( fT9l ) gives the velocity 
along the shell. 



v\ l9 — sin^cos 9 

V - — 

2 \ 1 - cos 61 



(A.4) 



On the other hand, the surface density of the layer is calculated as 
follows. First, we obtain the mass injection rate AM (into a solid 
angle defined by an impact parameter r) of the stellar winds. 



AM - rhxil - COS0) . 



EM, 



Anhn^clRl 



cos^9{l - cos 9y 



smd(9 - sin cos 0) 



(A.8) 



Let us now calculate the emission measure of the stellar 
winds. Since we assume the same flow parameters of both com- 
ponents of the binary system, then the emission measure of the 
winds is given by 

EM„ = 2 f ° nldz = 2 1-^]' f" , (A.9) 



and integrating one obtains 



EM„ = 



(r/Ro) 



n 

2 1 + (r/Roy 



+ arctg — 



(A. 10) 



(A. 5) Thus, the total emission measure is given by EM - EMs + EM„. 
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A.1. The emission measure EM as function of tine stagnation 
point radius Rq 

Let us now fix the impact parameter r, and investigate the be- 
havior of the emission measure EM as function of the stagnation 
point radius Rq. The total emission measure along a given line 
of sight can be written as 



EM : 



47T^ nf-clr' 



fs + U 



where 



and 

''-\{- 

4 \ vi 



co&^e(\ -coser 



sin 0(0 - sin cos 0) 



r/Ro ^ (Ro 

+ arctg — 

1 + (r//?o)2 ^\ r 



(A. 11) 



(A. 12) 



(A. 13) 



It is also possible to give the emission measure EMs 
(eq JA.SII ) in terms of r and Rq by substituting the trigonometric 
functions sin 9 = (r/T^o)/ a/I + {rl Rgf, c os 6* = 1 / ^j\ + (r/Ro)^, 
and - arctg(r//?o) into equation ( IA.12b . 

In Figure IATI we present fs and f^, as functions of r/Ro. 
We fix the impact parameter r and vary the stagnation point 
radius Rq. We note from the figure that, initially as r/Ro in- 
creases (which means closer binary systems), the emission mea- 
sure from the shell increases. While r/R^ is still growing, the 
emission measure reaches a maximum value (/, ^ 1 at r/Ro ^ 
1.3) and then it steadily decreases. In addition, it can be observed 
from the figure that /, » f^. at r/Ro ^ 1, and f <K /„, at very 
large (r/i^o 1) or very small (r/7?o » 1) stagnation point radii. 




0.01 0.10 1.00 10.00 100.00 1000.00 

(r/Ro) 



Fig. A.3. Integrated flux density at 43 GHz of the shocked layer 
over the nondimensional impact parameter r/Ro- In this exam- 
ple, we have assumed the parameters of the colliding winds of 
Model B2 (see Table 2). The behavior of the plot is described in 
the text. 



Consequently, with in this limit the optical depth of the shell 
does not depend on the value of r and scales as D^^ (being D the 
binary separation). 

On the other hand, for higher impact parameters [r/Ro » 1) 
fs^2/niRo/rf and 



Tt — Ay c ■ 



(A. 17) 



Thus, in this limit (t, oc r ^), the optical depth of the shell (for a 
given line of sight) scales as bc D^. 



A.2. The emission measure EM as function of the impact 
parameter r 

We now fix the stagnation point radius Rq and investigate the 
variation of the total emission measure with the impact parame- 
ter r. In this case, it is useful to write the total emission measure 
EM{^ EMs + EM„) as 



EM = 



An^m^ c^RI 



(A. 14) 



In Figure IA.2I we show the contribution of the shell EMs, 
the contribution of the winds EM^., and the total emission mea- 
sure EM as functions of the nondimensional impact parameter 
{r/Ro). We note from the figure that the emission measure is 
dominated by the shell {EMs ^ EM^.) at impact parameters 
(r/Ro) - 1. However, at very low {r/Ro «; 1) or very high 
{r/Ro s> 1) values of the impact parameter, the contribution of 
the stellar winds to the total emission measure is more important 
{EMs « EM,,). 

Given the emission measure EMs of the shocked layer, the 
optical depth Tj can be obtained by 



where Av = OTj;i'(v)/(4;7r^;r!^Cj). 



(A. 15) 



For lines of sight with impact parameters r/Ro <sc 1, it can 
be shown fliat ^ 3/8 (r//^o)^ and 



3 Ay 

" 8 R^ 



(A. 16) 



A.3. Flux density Sy as function of the binary separation D 

In this section, we investigate the dependence of the flux den- 
sity Sy with the binary separation D. Figure lA3] shows the flux 
density 5 43 GHz of the shocked layer as function of the nondi- 
mensional impact parameter r/Ro. We have adopted the wind 
parameters of Model B2 (see Table 2). It can be seen from the 
figure that at the beginning the flux grows as r^, which suggests 
that the emission comes from an optical thick disk. Eventually 
(at r/Ro - 40), however, the shell becomes optically thin and the 
flux tends to a constant value. The main contribution to the flux 
comes from the optically thick region of the shell, and therefore 
the emission from the optically thin region can be neglected. 

Let r,„ be the transition impact parameter at which Ts{r„) = 
1. Since EMs ^ EM„, for lines of sight with impact parame- 
ters r < r,„, we assume that the optical depth is dominated by 
the shell. (For simplicity, we are not considering those impact 
parameters r «; Rq for which EMs ^ EM^,; see Figure [ATI ) 
It follows that the optical depth T{r,„) ^ {2 / 7T)Ay{Rl/ rf„) ^ 1 
(eq. IIA.16II ). and then the size of the optically thick region in- 
creases as r„, oc R^^'. Finally, since the total flux is dominated 
by the emission from the optically thick region (r < r„,), one 
can assume that the flux Sy cc r^^. Consequently, we predict from 
our model that the emission from the shocked layer scales as 
Sy cc D^l^ . Detailed calculations (as in the example shown in 
Fig. IA.4I ). which consider both the contribution of the optically 
thin region to the emission from the shell and the contribution by 
the stellar winds to the total optical depth, give a small deviation 
from this prediction, that is, S'y oc D*'^'^. 
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Fig. A.4. Predicted flux density at S 43 ghz from a binary system 
as function of the binary separation D. We have adopted the stel- 
lar wind parameters of Model B2 (see Table 2). From this exam- 
ple, we found that the flux scales as Z)"-^^. 
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